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ABSTRACT 

We present the first results from the SWARMS survey, an ongoing project to identify compact 
white dwarf (WD) binaries in the spectroscopic catalog of the Sloan Digital Sky Survey. The first 
object identified by SWARMS, SDSS 1257+5428, is a single-Hned spectroscopic binary in a circular 
orbit with a period of 4.56 hr and a semiamplitude of 322.7 ± 6.3kms~^. From the spectrum and 
photometry, we estimate a WD mass of 0.92^0 32 1^©- Together with the orbital parameters of the 
binary, this implies that the unseen companion must be more massive than 1.62t°J0 M0, and is in 

all likelihood either a neutron star or a black hole. At an estimated distance of 481; i[g pc, this would 
be the closest known stellar remnant of a supernova explosion. 

Subject headings: binaries: close, spectroscopic — white dwarfs — supernovae: general 



1. INTRODUCTION 

The vast majority of stars with masses below 8M0 end 
their lives as white dwarfs (WDs). Left to their own de- 
vices, WDs will just cool and fade away, but WDs that 
form part of close binary systems can interact with their 
companions, steering away from the normal course of iso- 
lated stellar evolution and giving rise to a plethora of 
astrophysical phenomena. When the companion of the 
WD is another compact object (a second WD, a neutron 
star [NS], or a black hole [BH]), we refer to these systems 
as compact WD binaries, or CWDBs. The orbital evo- 
lution of detached CWDBs is governed by the ernission 
of gravitational waves (GWs, IPeters fc Mathewsl IT963I : 
lPaczvhsOll967f ). For systems with small enough sepa- 
rations, loss of angular momentum leads to a merger of 
the two components in a time that is short enough to 
be of astrophysical interest. This kind of CWDBs con- 
stitute the largest population of Galactic GW sources 
(jEvans et al.lll987h . domin ating the GW fo reground for 
future missions hke LISA (jNele man illml). They can 
also provide important constraints on the still poorly un- 
derstood process_ofconimon-envelope evolution in binary 
systems (jNelemans k, To ut 2005j). 

But perhaps the most interesting aspect of CWDBs 
is the outcome of the merging process at the end of 
their binary evolution. For systems composed of two CO 
WDs whose combined mass is above the Chandrasekhar 
limit, the merger might trigger a thermonuclear run- 
away a nd lead t o a T ype la supernova (SN) explosion 
(|Iben fc Tutu kOTlll984t [Wcbbink 1984). This possibihty, 
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often referred to as the double degenerate WD (DDWD) 
SN la progenitor scenario, is the only theoretical model 
that naturally explains the absence of H in the nebular 
spectra of Type la SN (see lLeonar dl [20(T7l and refer ences 
therein), and it has motivated a number of searches for 
suitable candidate systems. Th e most comprehensive o f 
these searches, the SPY survey (jNapiwotzki et al.ll2001h . 
examined the spectra of ~ 1000 WDs to look for radial 
velocity (RV) shifts in the characteristic absorption lines. 
The total number of detached DDWD s found to date by 
SPY and other surveys is around 100 (jNapiwotzki et al.l 
I2004D , but the periods and mass estimates have o nly been 
pubhshed for 24 systems (jNelemans et al.|[2005D . So far, 
none of these systems clearly fulfills the requisites to be 
a SN la progenitor. 

In this paper a nd in a companion publication 
(jMuUallv et al.|[2009[ ). we present the first results from 
SWARMS, the Sloan White dwArf Radial velocity data 
Mining Survey. The aim of SWARMS is to mine the 
WD catalog in the spectro scopic data base of the Sloan 
Digital Sky Survey (SDSS, j^rkeTaDllQOS) in search of 
CWDBs. Our ultimate goal is to find the DDWD pro- 
genitors of Type la SNe, or at least put constraints on the 
rate of DD mergers in the Galactic disk that can be com- 
pared with measurements of the local Type la SN rate. 
This will allow us to assess the viability of the DDWD 
progenitor scenario for SN la. This paper is organized 
as follows. Our data mining strategy is briefly described 
and compared to the SPY survey in Sectionj^l using as an 
example the CWDB SDSS J125733.63+542850.5 (hence- 
forth, SDSS 1257+5428), the first object discovered by 
SWARMS. In Section [31 we present the follow-up ob- 
servations of SDSS 1257+5428 that we have conducted 
in order to determine its orbital and spectral parame- 
ters. In Section |31 we discuss the implications of our 
findings for the nature of the unseen companion of SDSS 
1257+5428. We summarize our conclusions and outline 
the future prospects for our survey in Section [5l 

2. SWARMS: MINING THE SDSS DATA BASE FOR 
COMPACT WD BINARIES 

2.1. General Strategy 
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The firs t version of the S PSS W D catalog was pub- 
lished by iKleinman et al.l (120041). T h e mo st recent 
version, presented in lEisenstein "eFaLl pOOel ) (h ence- 
forth, E06), corresponds to Data Release 4 of the SDSS 
(jAdelman-McCarthv et ani2006D , and it is currently the 
largest catalog of its kind, with 9316 spectroscopically 
confirmed WDs. This number should be roughly doubled 
in the forthcoming DR7 version. Many publications have 
analyzed the SDSS spectra of WDs and related objects 
like subdwarf stars and cataclysmic variables (see E06 for 
a list of references), but so far all this work has focused on 
the final co-added spectra. All the SDSS spectra arc in 
fact taken in t hree or more exp osures, approximately 15 
minutes long (jStoughton et al.l [2002). in order to facili- 
tate cosmic ray rejection. The exact number of exposures 
and the time lapsed between them varies from plate to 
plate, depending on the observing conditions at the time, 
but exposures are usually taken back to back during th e 
same night. Starting with DR7 (|Abazaiian et al.|[2009[ ). 
the individual exposures are available separately, which 
opens the domain of time resolved spectroscopy in the 
SDSS data base. 

The goal of SWARMS is to take advantage of this sit- 
uation in the context of the WD spectra, looking for RV 
shifts among different exposures of the same object that 
may be the hallmark of a massive companion. To fa- 
cilitate the search, we have started by considering only 
WDs with 'dominant type' DA or DB in the E06 catalog, 
i.e., objects that have a relatively simple spectrum with 
absorption lines from either H or He. From this list, we 
have discarded any objects with excess flux in the red 
part of the spectrum, which are usually classified as de- 
tached WB+M star binaries (|Silvestri et al.ll2006[ ). This 
initial triage yields 7956 objects, which we have examined 
for RV shifts using several automated techniques supple- 
mented by visual examination. A detailed description of 
our search strategy and the implications for complete- 
ness, the derivation of merger rates, and other related 
issues will be the subject of a forthcoming publication. 

2.2. SWARMS vs. SPY 

Before going on to discuss specific results, it is in- 
structive to outline the fundam ental differences betwee n 
SWARMS and the SPY survey (jNapiwotzki et al.ll200ll ). 
At a very basic level, both surveys are complemen- 
tary, because the SPY objects were examined for RV 
shifts from the VLT telescopes in Paranal (Chile) and 
have therefore mostly southern declinations, while SDSS 
data arc taken from the Apache Point Observatory 
(APO) in New Mexico and have mostly northern de- 
clinations. The SPY survey achieved an excellent RV 
accuracy of ~ 2kms~^ by using the high-resolution 
UVES spectrograph (i? = 18500), but the demands of 
high-resolution spectroscopy imposed a stringent mag- 
nitude cutoff of B < 16.5 on the SPY WDs, even for 
the large collecting area (8.2 m) of the VLT Kueyen 
telescope (jNapiwotzki et al1l200lD . By comparison, the 
resolution that can be achieved by SWARMS with the 
spectrographs on the SDSS 2.5m telescope (i? = 1800, 
lYork et ani200Gf ) is much lower, around 170kms~^. But 
RV resolution is not the most important factor to find po- 
tential DDWD SN la progenitors. The most interesting 
systems are those with relatively massive components, 
or tight orbits, or both, which are expected to have RV 
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Fig. 1.— The three individual exposures of SDSS 1257-1-5428 
from the SDSS DR7 archive around the H/3, H7 and H(5 lines. 

shifts higher than ~ 150 km s"^ (jNapiwotzki et al.ll2004f) . 
Moreover, by foregoing the need for high-resolution spec- 
troscopy we can examine dimmer WDs and hence many 
more objects. Through a combination of centroid fits and 
cross-correlation techniques, we have been able to achieve 
an accuracy of ~ 120 km s~^, detecting RV shifts between 
consecu tive SDSS ex posures in a WD with a g magnitude 
of 18.9 (jMuUallv et a l. 2009). Thus, SWARMS can ex- 
amine roughly an order of magnitude more objects than 
SPY, albeit at a much lower resolution. We can summa- 
rize our comparison by saying that SPY is better suited 
for a systematic study of the properties of DDWD sys- 
tems, but SWARMS is better positioned to identify in- 
teresting objects, particularly the pre-mergers and SN la 
progenitors. 

2.3. The first SWARMS object: SDSS 1257+5428 

To illustrate the kind of objects that can be found by 
SWARMS, we present here results for SDSS 1257-f 5428, 
the first CWDB identified by the survey. In the E06 
catalog, SDSS 1257-^5428 was listed as a DA WD with 
a g magnitude of 16.8 (see Tabic [T]). There are three 
individual exposures of this WD in the DR7 SDSS data 
base. The WD is identified as a binary by the large shift 
in the centroid of the H lines between exposures and 

1 (taken on the night of October 3 2003) and exposure 

2 (taken on the following night; see Figure [1]). The shift 
has a value of 7.9 A at the line, corresponding to a 
RV difference of 487 km s"^ 

3. FOLLOWUP OBSERVATIONS OF SDSS 1257-1-5428 

3.1. Data acquisition and reduction 

After a WD in the SDSS catalog is identified by 
SWARMS as a member of a binary system, follow-up 
observations are necessary to determine its period and 
amplitude, and obtain high signal-to-noise spectra that 
can be used to measure its effective temperature T^ff 
and gravity logg, calculate its mass, and constrain the 
nature of its companion. In the case of SDSS 1257-1-5428, 
the follow-up observations were conducted with the Dual 
Imaging Spectrograph (DIS) at the 3.5m ARC telescope 
at APO. A total of 23 spectra were taken on the nights 
of February 5, 6, and 14, 2009. We used the B1200 grat- 
ing with a 1.5" slit to achieve a 0.62 A dispersion and a 
FWHM 1.8 A resolution. The integration time was 10 
minutes for all spectra, bracketed by exposures of He, Ne 
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TABLE 1 

Observational properties of SDSS 1257+5428 



Parameter 



Value 



Parameters from SDSS 



RA (J2000) 


12h 57m 33.63s 


Dec (J2000) 


+54° 28' 50.5" 


u (Mag) 


17.355 ± 0.009 


a (Mag) 


16.844 ± 0.004 » 


r (Mag) 


16.732 ± 0.005 


i (Mag) 


16.722 + 0.005 


z (Mag) 


16.762 + 0.010 


Spectral type 


DA b 


T.J J (K) 


8750 + 25 


log{g) 


9.0 + 0.005 1^ 



Orbital parameters (APO) 



Ka (kms-l) 


322.7 + 6.3 


P(hr) 


4.5550 + 0.0007 


ro(JD) 


2454868.1182 + 0.0006 


7^ (kms-i) 


-28.9 + 4.6 




44.4 



Spectral and Photometric Parameters (APO, SDSS) 



Teff (K) 


~ 9000 <i 


log{g) 


~ 8.5 


Ma{Mq) 


f\ Qo"t"0-28 c 
"■='^-0.32 


D (pc) 


48+10 o 


tCooi (Gyr) 


2.0 + 1.0 ° 



Companion Parameters 



MBSin{i) (Mq) 


1 62+"-^" 


Mb for i = 60° (M©) 


2 in+" " 

^- J^^-0.30 


tMerge for i = 90° (Myr) 


.,-,+342 


tMerge for i = 60° (Myr) 


26712.155 



DR7 values. 

Listed autof it results from E06. 
<^ For 23 data points (19 d.o.f). 

Spectral parameters from the APO fits arc listed without 
toleran ce ra nges due to the low quality of the fits, see 
Section 13.31 

See Section [3.3l for an explanation of the tolerance ranges 
on Ma and D. 

and Ar arc lamps to ensure a proper wavelength calibra- 
tion. Flux was calibrated each night by taking several 
exposures of the spectrophotometric standard Feige 67 
(|Okdll990n . The data were reduced using standard long- 
slit IRAF routines (|Todvlll993n . A prominent Hg emis- 
sion line at 4358 A, courtesy of the inhabitants of White 
Sands, NM, confirmed that the wavelengths are cali- 
brated to better than the instrumental resolution: the 
standard deviation of the line over 8 spectra (one night) 
is 0.1 A. The final spectra span between 3770 and 5030 
A, and clearly show all the lines in the Balmer series from 
H/3 to HiQ. 

3.2. Orbit 

We measured the RV shift of each APO spectrum by 
fitting the centroid of the line after applying the stan- 
dard Solar System barycentric correct ions. We used a 
Loren tzian profile with a Gaussian core ijThompson et all 
I2004D . adjusting the fit parameters with MPFIT, a ver- 
sion of the Leve nberg-Marquardt algorithm adapted for 
IDL (|Markwardt,i2009il . The resulting RV curve is well fit 
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Fig. 2.— RV curve for SDSS 1257+5428, folded in phase with 
the best-fit period. 
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Fig. 3. — Goodness-of-fit parameter (x'^) as a function of period 
for the RV curve of SDSS 1257+5428. 

by a circular orbit of the form Ka sin((27r(t — Tq) / P)) + 
7^, where P is the period, Ka the semiamplitude, Tq a 
temporal offset marking the zero-point of the curve, and 
7^ a constant RV offset (see Figure[2]). We found that the 
nine day baseline of the APO spectra strongly constrains 
the period of the binary to P = 4.5550 ± 0.0007 hr, (Fig- 
ure[3|), with a semiamplitude of Ka = 322.7±6.3kms^^ 
(see Table[l]for a complete list of the orbital parameters). 

3.3. Spectrum 
3.3.1. Spectral fits 

The autof it results from E06 give Te// = 8750 ±25 K 
and \ogg = 9.0 ± 0.005 for SDSS 1257+5428, but these 
parameters are derived from the co-added DR4 SDSS 
spectrum of an object that was not known to be a short- 
period binary. In order to make the best possible de- 
termination of Teff and log g, we have co-added the 23 
APO spectra after removing the RV shifts. The result- 
ing spectrum has a S/N of 160, much higher than the 
co-added SDSS spectrum (S/N=30). 

We have used the MPFIT routine to fit the co-added 
APO spectrum with the latest (June 200 8) version of 
the K oester atmosphere models descrbed in lFinlev et al.l 
(|1997| ). adjusted with a first order polynomial. This 
spectral fitting has proved a challenging endeavor. The 
raw (unreduced) parameter has two distinct minima 
in the T^ff, logg plane: one around T^ff = 35000 K, 
log g = 7.5 and one around T^/f = 8500 K, logg = 8.5 
(see Figure lU left panel). The hot solutions provide 
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Fig. 4. — Constraints on T^f f and logg for SDSS 1257+5428. Left: surface contours of the raw parameter obtained from fitting 
the co-added APO spectrum of SDSS 1257+5428 to the most recent version of the models by IFinlev et al.l 1)19971 ). The best-fit solution 
{Tf^ff = 9000 K, logg = 8.5, see text and Figure [Sj is indicated by a large red circle. The models used to derive the confidence ranges 
for the WD mass are indicated by smaller blue circles, joined by straight lines. Right: co mparison between the DR7 SDSS photometry 
of SDSS 1257+5428 (red circle) and the photometric models of lHolberg fc BergeronI 1)2009 ). The cooling tracks are labeled by log 3, with 
points at the same Tf,ff represented by dotted lines. 

slightly lower values than the cold solutions, but nei- 
ther fits are very impressive (see Figure [5] for the best- 
fit cold solution - the visual quality of the best-fit hot 
solution is similar). This degeneracy between cold and 
hot spectral solutions is common in DA WDs, and it 
can usu ally be broken by inc luding photometric infor- 
mation (|Kleinman et al.ll2004[ ). In the right panel of Fig- 
ure H we compare the DR7 SDSS photometry of SDSS 
1257+5428 (Table [B with the photometric models of 
IHolberg fc Bergeronr(|2006n . The hot solutions can be 
discarded with a high degree of confidence. The cold 
solutions come very close to reproducing the SDSS pho- 
tometry, with the models along the cold minimum 
being offset from the data by only ~ —0.15 magnitudes 
in u — g. This kind of discrepancy is not uncommon in 
comparisons betwee n the SDSS photometry of D A WDs 
and cooling models (jHolberg fc Berg eron 2001 see also 
Section 13.3.21 below) . Includin g the full reddening cor- 
rection of iSchlegel et al.l (|1998f ) changes the u — g color 
by —0.025 magnitudes and the g — r color by —0.018 
magnitudes, which reduces the discrepancy by a small 
amount. 

Even though it is not exactly at the bottom of the lo- 
cal minimum, we have found that the spectral model 
with Teff = 9000 K and logg = 8.5 (Figure O provides 
the best match to the continuum and the cores of the 
high order Balmer lin es, which sh ould be more sensitive 
to the value of logg (jFinlev et a l. 1997). We consider 
this model our 'best fit', although it cannot be said that 
the spectral adjustment shown in Figure [5] is particu- 
larly good. Nevertheless, the values of T^/f and logg 
that we obtain are reasonably close to those found by 
E06 using the co-added SDSS spectrum. This agreement 
between different analysis techniques applied to different 
data sets of the same object is encouraging, and supports 
our choice of spectral model. 



3.3.2. Why are the fits so poor? 

The failure of theoretical models to reproduce the 
details of the spectrum and photometry of SDSS 
1257-1-5428 is not entirely unexpected. Even the most 
updated versions of the Koester atmosphere models have 
problems for WDs with T^ff < 12000 K. At these tem- 



peratures, the contribution of neutral and charged per- 
turbers to the equation of state is of the same order of 
magnitude, and the size of the c qnvective enve lope of 
the WD has an abrupt increase (jKoester et al.l [2008) . 
These effects have been put forward as possible causes 
for the observed increase of the aver a.ge effective gravity 
for WDs in this tem perature range (jKepler et al.l [20071 : 
[Bergeron et al.|[2007|) . The discrepancies shown in Figure 
[5l however, are rather large, and they should not be writ- 
ten off as shortcomings of the models before considering 
alternative explanations. 

The high S/N of our co-added APO spectrum reveals 
some unusual features in the continuum around 4600 A 
and 4750 A. The features are not present in the pho- 
tometric standard Feige 67, which indicates that the 
flux calibration is not at fault. None of the calibration 
frames, flat fields or biases have fiaws that could have 
caused these features. Other WDs that w e observed with 
the s ame spectrograph in the same mode (jMuUallv et al.l 
I2009D do not show them, but these are fainter objects 
with worse S/N in the co-added spectra. At this stage, 
we have no explanation for the presence of these localized 
features, but they are clearly not large enough to infiu- 
ence the mismatch between the models and the observed 
spectrum. 

It is possible that the spectrum of SDSS 1257+5428 is 
simply distorted by optical emission from its binary com- 
panion. The co-added SDSS spectrum shows no traces 
of line emission or a substantial excess in the red spec- 
trograph, clear signs of the presence of a nondegenerate 
companion that would have prompted a DA-I-M classifi- 
cation for SDSS 1257-^5428 by E06 (this is ruled out by 
the RV curve as well, see the discussion in Section [4]) . 
A WD companion is more difficult to discard from the 
point of view of the spectrum alone. We have examined 
the individual APO spectra that are closer to quadrature, 
when the relative RV shift of the individual components 
should be larger, but found no sign of artifacts or asym- 
metries in the Balmer lines that might indicate the pres- 
ence of a second DA WD of similar temperature. A sub- 
stantially hotter WD of any spectral class would rapidly 
overwhelm the spectrum of SDSS 1257+5428, having a 
large impact on the SDSS photometry shown in Figure 
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Fig. 5.— Best-fit spectral solution for SDSS 1257+5428: T^// = 9000 K, logg = 8.5. 
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Fig. 6.— Comparison between the SED of SDSS 1257+5428 from 
the five -band DR7 SDSS photomet ry (red circles) with the predic- 
tions of lHolberg &: BereeronI 120061) for WD models along the cold 
minimum, T^j j between 8000 and 9000 K and logg between 
8.0 and 9.0 (shaded area). 



131 The only possibility left is a non-DA WD of similar 
or lower temperature with a rather featureless spectrum. 
In Figureiniwe display the spectral energy density (SED) 
of SDSS 1257 +5428 together with the ra nge covered by 
the models of lHolberg fc BergeronI (pOOl ). The SEDs of 
models along the cold minimum are roughly consis- 
tent with the five-band photometry of SDSS 1257+5428. 
There might be a hint of an excess in the z band, but in- 
frared data would be required to confirm this possibility. 
On the other hand, all the models seem to overprcdict 
the flux in the u band, an effect that was already ap- 
parent in the color-color plot shown in Figure [31 Since 
this is a flux deficit instead of an excess, it is hard to ex- 
plain by invoking th e prese nce of a companion. In fact, 
IHolberg fc BergeronI (|2006D find similar problems in the 
u band for other cold DA WDs (see their Figure 6). 

We conclude that the shortcomings of the WD models 
in this temperature range are indeed the most likely ex- 
planation for the poor quality of the spectral fit shown 
in Figure O A cold WD companion cannot be ruled out 
based on the spectral information alone, but neither is 
there any evidence in the observations that requires the 
presence of such a companion. We shall revisit this pos- 
sibility in Section [31 

3.3.3. Estimating the mass of the primary WD 



The issues with the models that we have discussed re- 
quire that we be extremely careful with the interpreta- 
tion of the spectral fits. In particular, we cannot justify 
the calculation of confidence ranges around the fitted pa- 
rameters in the usual way (i.e., based on the statis- 
tic). But we are not interested in the confidence ranges 
around T^ff and log 5 per se. Rather, our goal is to use 
the spectrum of SDSS 1257+5428 to measure its mass 
and constrain the nature of its unseen companion. 

An estimate for the mass of SDSS 1257+5428 can be 
obtained by examining the performance of the spectral 
models around t he cold minimum in x^ space. The 
cooling cu rves of [Fontaine et al.l (120011) included in the 
models of IHolberg fc BergeronI (|2006D give a value of 
Ma = 0.92 M0 for our best-fit spectral solution, which 
is probably a good initial guess for the mass of SDSS 
1257+5428. Along the cold x^ minimum, the spectral 
solutions with Te// = 8250, 8500, and 8750 K; and 
log g = 8.0, 8.5, and 9.0 (indicated by the small blue cir- 
cles in the left panel of Figure [3]) appear to bracket the 
effective gravity values that are in best agreement with 
the higher order Balmer lines in the spectrum, as shown 
in Figure [71 Any models with log g outside this range 
result in extremely poor fits and must be discarded, in- 
cluding the models with lower T^f / and lower log g that 
seem to be preferred by the photometry (right panel of 
Figure [3]) • Since we have already ruled out a DA WD 
companion of similar temperature, we expect the shape 
of the higher order Balmer lines to be a reliable probe 
of the mass of SDSS 1257+5428, even if the underlying 
continuum is somewhat distorted. 

The spectral models shown in Figure [3 translate into 
masses between 0.60 and 1.20 M0, which we will adopt 
as conservative lower and upper limits to Ma- It is un- 
fortunate that the poor performance of the spectral mod- 
els in this temperature range does not allow for a more 
accurate measurement of Ma, but we believe that our 
estimate is reasonable, and probably the best that can 
be done with current tools. We note that the models 
we contemplate here span a range of log g values around 
the best-fit solution (±0.5) that is comfortably larger 
than the reported increase for Teff < 12000 K (~ 0.2, 
[Kepler et aLll2007l : [Bergeron et al.ll2007[ ). Also, our lower 
limit to Ma is close to the mean mass of DA W Ds with 
Teff > 12000 K (0.58 M^. [Kepler et all [2001 . which 
would be the best guess for the mass of a WD in ab- 
sence of any spectral or photometric information. For 
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Fig. 7. — Comparison between the high order Balmcr hnes (H^, 
He, and Hg) in the spectrum of SDSS 1257+5428 and the three 
models along the cold minimum indicated with small blue circles 
in the left panel of Figure |4] Although higher order lines (Hg 
and Hio) arc present in the spectrum, the S/N is too low for a 
meaningful comparison with the models. 

all these reasons, wc are confident that the true mass of 
SDSS 1257+5428 lies within the range of values that we 
propose here. For comparison purposes, the best-fit hot 
spectral model {Tejf = 35000 K, logg ~ 7.5) translates 
into a value of 0.45 M0 for Ma- With the set of limiting 
spectral models around the best-fit c old solution shown 
in Figure HI the iHolberg fc BergeronI curves yield a cool- 
ing age of 2.0 ± 1.0 Gyr. The absolute g magnitude is 
13.4^0 4, which results in a distance of D = '^^-w PC for 
SDSS 1257-1-5428. 

4. THE COMPANION OF SDSS 1257+5428: NEUTRON 
STAR OR BLACK HOLE? 

Combining our estimate for the mass of SDSS 
1257+5428 (Ma = 0.921^:^^ Mg) with the values of P 
(4.5550 + 0.0007 hr) and Ka (322.7 + 6.3 km s"!) derived 
in Section ESI we obtain MBsin(i) = 1.62tSj;25M0 for 
its unseen companion, with Mb and the inclination angle 
i being degenerate as in all single-lined binaries. A plot 
of the companion mass as a function of cos(i) is shown in 
Figure [H We stress that the properties of the RV curve 
alone require a massive companion, regardless of the es- 
timated value for Ma- Mb must be more massive than 
0.66 Mq even for Ma = (i.e., assuming the value of Ma 
is negligible for the dynamics of the system), and more 
massive than 1.08 Mg for Ma = 0.3 M©. With any rea- 
sonable range of masses for the primary, a nondegenerate 
stellar companion would have a spectral type G or ear- 
lier, which is clearly incompatible with the observations 
of SDSS 1257+5428. The companion must therefore be 
a compact object, either another WD, a NS or a BH. 

Could the companion be a WD? We have not been 
able to discard this possibility from the point of view of 
the optical spectrum alone, but the properties of the RV 
curve make it very unlikely. The companion could only 
be a WD if our mass estimate for Ma was completely 
off, and even then, it would have to be a very massive 
object (above 1.0 M© for any reasonable value of Ma), 
cool object, without any prominent absorption lines (see 
discussion in Section r3.3.2p . Such massive WDs are ex- 
tremely rare. But we have found no evidence to indicate 
that the shape of the Balmcr lines in the spectrum of 
SDSS 1257+5428 is distorted in any way, so there is no 




cos(i) 



Fig. 8.— Mass function for the companion of SDSS 1257+5428. 
The solid blue curve represents the values obtained using the best 
estimate for the WD mass {Ma = 0.92 M©), and the striped region 
is bound by the upper and lower limits to the WD mass obtained 
from the spectrum (Ma = 0.60 and 1.20 Mq). For comparison, 
the mass functions for Ma = 0.3 Mq and Ma = are also shown 
(dash-dot ted and dashed red plots). T he range of masses for known 
NSs from lLattimer fc PrakashI II2007I ) has been highlighted with a 
shaded region, taking into account the recent re-evaluation of the 
mass of PSR 0751+1807 bv lNice et all (120081) . 

reason to doubt our conservative estimate for Ma- If this 
estimate holds, assuming the smallest possible value of 
Ma (0.6 M0), the orbit would have to be nearly edge-on 
(z > 82°) for Mb to be below the Chandrasekhar limit 
Mch- This inclination is expected to happen randomly 
in only 14% of binary systems, but SWARMS does have a 
strong observational bias towards finding high inclination 
systems. Even then, the companion would be more mas- 
sive than the largest known WDs (1.3 3 to 1.35 M0, see 
iKepler et ani2007i : [Barstow et al.lll995f) . In this unlikely 
circumstance, the SDSS 1257+5428 system would have a 
combined mass at least 43% larger than Mch, but might 
not be a SN la progenitor if the com panion turns out to 
be an ONe WD instead of a CO WD (|Garcia-Berro et all 
Il997f l. Considering all the evidence, the most likely pos- 
sibility by far is that the companion of SDSS 1257+5428 
be the stellar remnant of a supernova explosion, either a 
NS or a BH. We note that our estimated distance range 
places this objec t closer to the Sola r System than any 
other known NS (jPosselt et al.ll2007f ). 

The puzzle of the nature of the companion of SDSS 
1257+5428 cannot be solved with the observations that 
we present in this paper, and must be the subject of fu- 
ture work. It is, however, interesting to speculate about 
the possibilities. For our best-fit value of Ma (0.92 M©), 
Mb is below 1.86 M© (the largest masured mass fo r a 
NS, see iLattimer fc PrakashI [20071 : iNice et ani2008f ) for 
i > 67°, which has a random likelihood of 39%. In this 
case, the mass of the WD and the circularity of the orbit 
would place the system in the class of 'intermed iate-mass' 
WD+NS binaries (see Table 1 in lStairsl [2004l). Current 
models for binary stellar evolution predict that these sys- 
tems undergo un stable mass transfer during a common- 
envelope phase (jStairsI l2004( ). which results in the NS 
becoming a mildly recycled millisecond pulsar (MSP). 
This scenario would explain the high inferred mass for 
the companion of SDSS 1257+5428, well above the aver- 
age for a NS, as a byproduc t of the mass transfer process 
(jLattimer fc Prakashl[2007l ). If the NS is indeed an MSP, 
we expect the magnetic field to be low (~ 10^ Gauss), 
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the pulsar lifetime to be large (~ 1 Gyr), and the open- 
ing angle to be wide (poss ibly many tens of degrees, see 
iPhinnev fc Kulkarnilll99^ . Under these conditions, the 
prospects for detecting such a nearby MSP are good, 
although we note that the cooling time for the WD is 
also of the order of Gyr (see Section [X^ . and the MSP 
might have lost a large part of its magnetic field. For 
Ma = 0.92 Mq and inclination angles below 67° (61% 
random likelihood), the companion would probably be a 
stellar mass BH. In this case, we would not expect any 
kind of direct emission from it. 

No radio or X-ray source appears at this location in 
any of the indexed astronomical catalogs, but this does 
not preclude the existence of a faint counterpart to the 
companion of SDSS 1257+5428. This part of the sky 
has only been shallowly surveyed for MSPs in the ra- 
dio. The most stringent limits are probably from the 
Green Bank 140 ft 350 MH z survey, which h ad a nomi- 
nal sensitivity of 12-15 mJy (jSaver et al.lll997[) . In the X- 
rays, this location has never been observed with Chandra 
or XMM-Newton, and a faint nearby source could have 
easily escaped detection by the ROSAT All-Sky Survey. 
A systematic cross-correlation between stellar sources in 
SDSS (including all the WDs from E06) and the R OSAT 
catalogues was performed bv lAgiieros et al.] (|2009D . who 
found no counterpart to SDSS 1257+5428. 

One interesting implication of the companion of SDSS 
1257+5428 being a NS or BH is that the system should 
have received some kind of kick from the SN explosion. 
In principle, the spatial velocity of SDSS 1257+5428 
can be determined by the temporal average of the non- 
gravitational Doppler shifts in the spectrum, which gives 
the radial component, and the proper motion measured 
by SDSS {fi = 0.049 mas yr~^), which gives the compo- 
nent on the plane of the sky. Taking our distance es- 
timate, the proper motion translates into a transverse 
velocity of llj^4kms~^. The radial component is more 
difficult to estimate, because the constant component to 
the RV curve (7^ = —28.9 ± 4.6kms~^ in our fit, see 
Section 13. 2p is the combination of the true radial veloc- 
ity and the gravitational redshift of the WD. The value 
of the gravitational redshift depends on the WD mass, 
which we cannot measure with accuracy (see Section r3.3p . 
and the WD radius, which is also model dependent. A 
detailed estimate of the gravitational redshift for SDSS 
1257+5428 is outside the scope of this work, but for a 
massive ^ 1 Mp) WD, we expect it to be of the order of 
90kms-i (jWegner fc Reidlll991l) . which would require a 
radial velocity around —120 km s~^. The total spatial ve- 
locity would then be ~ 120 km s^^, mostly in the radial 
direction, which is compar able to the measur ed kicks for 
pulsars in binary systems (|Wang et al.ll2006[ ). 

Regardless of what the nature of the companion to 
SDSS 1257+5428 turns out to be, the system is clearly 
very interesting from an astrophysical point of view, and 
new observations should yield exciting results in the near 
future. From the measured orbital parameters, the sep- 
aration of the components must be small, with values of 
the scmimajor axis larger than 0.02 AU only for i < 17° 
(0.0093 AU at i = 60°). If pulsations from a NS compan- 
ion were detected, this could allow for a significant mea- 
surement of the Shapiro delay, as in PSR J0437— 4715 
(jvan Straten et aLll200in . The merging time t]\{erge of 



the SDSS 1257+5428 system due to GW radiation is 
< 511^^4^ Myr , with the uncertainty in this upper limit 
arising from the uncertainty in the determination of Ma- 
For the canonical inclination angle i = 60°, t Merge be- 
comes 267trf Myr. We note that only three WD+NS 
binaries with tMe.me < tnnh hie were known previously 
(IStairsll2004l : iKim et al.ll200l . and SDSS 1257+5428 has 
a shorter period and tMerge than any of them. The dis- 
covery of a fourth object in this class should prompt a 
revision of the estimated WD+NS merger rate s and their 
expec ted contribution to the GW background ()Kim et al.l 
[2OOI . 

The ultimate f ate of SDSS 1257+5428 is unclear. 
IKing et all (|2007[ ) have proposed the merging of mas- 
sive WDs onto NSs as a scenario for the origin of 
long-duration gamma-ray bursts without an accom- 
panying SN like GRB 060614 (|Gal-Yam et all I2OO6I : 
iDella Valle et al.] I2OO60 . but theoretical simulations for 
the final accretion phase in this kind of events have not 
been performed yet. The coalescence of a BH with a mas- 
sive W D may also lead to a gamma-ray burst (jFrver et al.l 
II999D , but the observational signature depends again on 
the details of the final accretion phase, which are poorly 
understood. More exotic astrophysical transients, such 
as may be observed by pres ent fPTF. [Ran et al.l [20091 ) 
and forthcoming (e.g. LSST. ITvezic et al.ll2008[ ) synoptic 
surveys, cannot be discarded as outcomes. 

5. CONCLUSIONS 

In this paper a nd a companion publication 
(jMuUallv et al.l l2009f ) we have presented the first 
results from the SWARMS survey, an ongoing project 
aimed at discovering and characterizing CWDBs in 
general and DDWD SN la progenitors in particular 
among the WDs in the spectroscopic SDSS data base. 
Further results and a complete description of the survey 
will be the subject of forthcoming publications. Our 
ultimate goal is to estimate a rate for the merger of 
DDWD SN la progenitors in the Galaxy, and to use 
that rate, in combination with the results from SPY 
and other surveys, to assess the viability of the DDWD 
progenitor scenario for Type la SNe. 

To demonstrate the capabilities of SWARMS, we have 
focused on SDSS 1257+5428, the first object found by the 
survey. This WD was identified as a member of a binary 
system from the SDSS spectra, and follow-up observa- 
tions were performed to measure its orbital and spectral 
parameters. The RV curve has a period of 4.5550+0.0007 
hr, with a semiamplitude of 322.7 ± 6.3 km s^^. The pri- 
mary WD is cold (~ 9000 K), which makes the spectral 
analysis somewhat challenging, but we have derived a 
conservative estimate of 0.92;to'32 for its mass. This 
implies that the unseen companion must be larger than 
1.62to 25M0, and is probably a NS or a BH. At a dis- 
tance of D = 48^}g pc, this would be the closest remnant 
of a SN explosion to the Solar System. We have discussed 
the implications of our discovery, and suggested future 
avenues of research on this object. This kind of massive, 
close binary is a perfect example of the objects that can 
be found by SWARMS and oth er similar surveys like 
MUCHFUSS (iTiUich et al l 120091) . We anticipate more 
such discoveries, which we will report in the literature. 
We also plan to set up an on-line database of SWARMS 
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objects. 
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